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NUMERICAL COMPUTATION OF A CERTAIN DIRICHLET 
SERIES ATTACHED TO SIEGEL MODULAR FORMS OF 

DEGREE TWO 

NATHAN C. RYAN, NILS-PETER SKORUPPA, FREDRIK STROMBERG 

Abstract. The Rankin convolution type Dirichlet series Df^g(s) of Siegcl 
modular forms F and G of degree two, which was introduced by Kohnen and 
the second author, is computed numerically for various F and G. In particular, 
we prove that the series Dfq(s), which share the same functional equation and 
analytic behavior with the spinor L-functions of eigenforms of the same weight 
are not linear combinations of those. In order to conduct these experiments a 
numerical method to compute the Petersson scalar products of Jacobi Forms 
is developed and discussed in detail. 



1. Introduction and Statement of Results 

For cusp forms F and G in the space Sk of Siegel cusp forms of degree 2 and even 
weight k on the full Siegel modular group, let Dp t c{s) be the Rankin convolution 
type Dirichlet series which was introduced in [9]. It was proved [9, first theorem] 
that Df^g(s) shares the same analytic properties and the same functional equation 
as the spinor L-function L(F, s) associated to a Hccke eigenform F in Sk- In fact, if 
G is in the spezialschaar then one has [9, second theorem] Df,g{s) = (Fi\Gi)L(F, s) 
(see Section 2 for notations). However, it was left open whether Df.g(s), for 
interesting Hecke eigenforms 1 F, G, is related to the spinor L-functions of the forms 
in Sk- In any case it cannot in general simply satisfy the same identity as before 
since, for example Dp^is), for F not in the Maass spezialschaar, has a simple 
pole at s — k with residue equal to the Petersson norm of F, whereas L(F, s) is 
holomorphic. In fact, a spinor L-function of a Hecke eigenform in Sk has a pole if 
and only if F belongs to the Maass spezialschaar [5] . 

The first interesting Hecke eigenform is in weight 20 and the space of such forms 
is 1-dimensional in weights 20 and 22, is 2-dimensional in weights 24 and 26 (which, 
remarkably, have rational Fourier coefficients and eigenvalues), is 3-dimensional in 
weight 28 and 4-dimensional in weight 30. The interesting eigenforms in weights 20 
up to 26 have rational Fourier coefficients (and eigenvalues) , whereas the eigenforms 
in weight 28 and 30 are conjugate under Galois conjugation, respectively. The aim 
of this article is twofold. We report the results of numerical experiments which we 
performed to study the series Dp,G f° r interesting eigenforms F and G in weights 
20 up to 30. Secondly, we develop and describe the necessary theoretical tools for 
our experiments. Necessary to the computations which we report about out in this 
note is also the ability to compute Siegel modular forms, their Hecke eigenvalues 
and their spinor L-series. We do this via a Sage [12] package written by the first two 
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authors [10] (with the aid of several people) and based on the algorithms developped 
in [15]. The code and the data which we used for our computations can be found 
under [11]. The interested reader may use these data to reproduce our experiments. 
Our computational results are summarized as follows. 

Main Theorem. For k = 20,22,24,26,28,30, let F and G denote simultaneous 
eigenforms of weight k that are cusp forms but not in the Maass spezialschaar. Then 
Df,g( s ) is n °t a linear combination of the spinor L-functions L(H,s) where H 
ranges over all simultaneous eigenforms (including non cusp forms) of weight k. 

This result is surprising since, as already pointed out, the series Df,g(s) shares 
the same functional equation and T-factors as L(F,s) and L(G,s) yet seems not 
to be related to the usual spinor L-functions. Note that the theorem is proved by 
numerical computations. However, we made careful error estimates to ensure the 
validity of our result. 

In fact, we actually prove more. Namely, we show in section 4 that all the series 
Df,g(s) and L(H, s), where {F, G} ranges over all subsets with one or two elements 
of the interesting eigenforms of weight k, and where H ranges over all eigenforms 
of weight k, are linearly independent. 

Investigating the (numerically approximated) coefficients of the series Df,g(s) 
more closely we find: 

Supplement to the Main Theorem. In the notations of the Main Theorem the 
coefficients of Df.g{ s ) ar e n °t multiplicative; that is, the (normalized) coefficients 
a(n) do not satisfy a(m)a(n) — a(mn) for all coprime pairs of integers m and n. 

We remark that this supplement implies that the series Df,g( s ), f° r F and G 
not is the Maass spezialschaar, is not an L-function in the sense of Sclbcrg, i.e., 
does not belong to the so-called Selberg class. Note also that this is the second kind 
of Dirichlet series which is attached to a Siegel modular form, has a meromorphic 
continuation, satisfies a functional equation but is not an L-function in the sense of 
Selberg, the other being the Kocher-Maass series. It would be interesting to have 
a deeper understanding of how these series Df,g(s) fit into the general framework 
of automorphic forms and their associated L-functions. 

In Section 2 we summarize the theoretical framework for our experiments. In 
Section 3 we develop the necessary tools for our computations and discuss com- 
putational issues. In particular, in order to obtain our numerical results we have 
to develop methods to compute Poincare series and Petersson scalar products of 
Jacobi forms. These methods are of independent interest; they are described in 
detail in Sections 3.1 to 3.3. In Section 4 we explain the proof of the main theorem. 
Finally, in section 5 we summarize our numerical results. 

2. Background and Notation 

The Siegel upper half space of degree d is given by the set 

U d := {Z = X + iY e C dx d : Z t = Z, and Y positive definite}. 

A Siegel modular form F : M d — > C of degree d > 1 is a holomorphic function 
invariant under the action of Sp(d, Z) on H d . In this paper we specialize to degree 2 
and so we have a Fourier expansion of the form 

F(12)= Y, a F (n,r,m)q n Cq' m , 

r 2 — 4mn<0 
n,m>0 

where q = e 2mT , q' = e 2mT ' and C = e 2mz (t,t' gi'.zeC). Such an F is a cusp 
form if the ap (n, r, m) = unless the quadratic form [n, r, m] = nx 2 + rxy + my 2 
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is strictly positive definite. We denote by Mk and Sk the space of Siegel modular 
forms and cusp forms of weight k and degree 2, respectively. 

For each integer I > 1 there is a Hecke operator T(£), a linear map that maps 
Sk to Sk (see [15, p. 386] for explicit formulas). We call a cusp form an eigenform 
if it is a common eigenform for all Hecke operators. To eigenforms we attach the 
spinor L-function as follows. If T(£)F = Xf^F, then define 



oo * 

(1) L(F,s) :=C(2s-2fc + 4)^ 



F,l 



This function has a meromorphic continuation to the whole complex plane with at 
most a pole, which is then at s = k and of order one. It satisfies the functional 
equation 

(2) L*(F,s) := (2n)- 2s T(s)T(s - k + 2)L(F,s) = (-l) k L*(F, 2k - 2 - s). 

In some of our computations we renormalize (2) by s n- s + k — | so that the 
corresponding new functional equation expresses an invariance under s 4 1 — s 
instead of s n- 2k — 2 — s. 

We use Sk,N for the space of Jacobi cusp forms of weight k and index N in the 
sense of [4, Chapter I]. Any <j) in Sk,N has a Fourier development of the form 

^( V )=^C^,r),T( r . 

Here the sum is over all N -admissible pairs, i.e. over all pairs of integers (D, r) such 
that D < and D = r 2 mod 47V. Recall that C<p(D,r) depends on r only modulo 
2N and that is a cusp form if C^(D, r) — unless D < 

As mentioned above, for an integer k we denote by Sk the space of Siegel cusp 
forms of weight k on the full Siegel modular group, and we use S^. and S ? k for, 
respectively, the subspace of Maass cusp forms and its complement. Recall that 
SI and Si are invariant under all Hecke operators. The first case of nontrivial 
Sj, occurs for weight k = 20, where Sl is one-dimensional. In Table 1 we list, for 
20 < k < 30, the dimensions of various subspaces of the full space of Siegel modular 
forms of degree 2 and weight k. 

Every element F in Sk admits a Fourier- Jacobi development 



2-kiNt' 

N>1 



with elements Fn of Sk,N, respectively. 

If and G are both elements of Sk we consider the Dirichlet series 

(3) ^, G ( S ):=C(2 S -2fc + 4)E^M, 

N>1 

which was introduced in [9] and which may be viewed as an analogue of the Rankin 
convolution for elliptic modular forms. Here (4>\ip), for Jacobi cusp forms <j) and ip 
in Sk,N, denotes their Petersson scalar product. More precisely, we have 

((j)\ip) = / </>(t, z)ip(r, z) v k ~ 3 e~ 4nNy2/v du dv dx dy (r = u + iv,z = x + iy), 
Jv 

where V denotes any fundamental domain for the Jacobi group SL(2,Z) J acting 
on H x C. 
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3. Method of computation 

3.1. Petersson scalar products. In order to compute numerically the coefficients 
of the L-series Dp,G we need a method for computing the scalar product of two 
Jacobi forms of the same weight k and index N . Since F and G are eigenforms, 
their Fourier coefficients (and hence also the Fourier coefficients of the Jacobi forms 
in their Fourier- Jacobi development) are algebraic numbers. The package [10] im- 
plements methods developed in [15] to provide these Fourier coefficients in an exact 
representation. It is thus reasonable to look for a method to express the Petersson 
scalar product of two Jacobi forms in terms of their Fourier coefficients. We propose 
the following method. 

For an iV-admissible pair a = (D,r) denote by P a the a-th Poincare series of 
Sk,N- This is the unique Jacobi cusp form which satisfies 

(<P\P a ) = C7 (a) 

for all forms </> in Sk,N- It is clear that the P a with a running through all admissible 
pairs span the space of cusp forms. Let n be the dimension of Sk,N, and let 
h = {h(i)} (1 < i < n) be a sequence of n admissible pairs. We use for the 
Gram matrix of the Ph(i) with respect to the Petersson scalar product, i.e. 

(4) T h ={C PHi) {h{j))). 

Note that the determinant of Th is different from if and only if the Ph(i) form a 
basis of Sk,N- For a form <f> in Sk,N, let 

lh {<i>) = {C rt> {h{j)),...,C f {h{n))) t . 

We then have the following formula for the scalar product of F and G. 

Theorem 1. Assume that the Poincare series Ph(i) A — * — n ) f orm a basis of 
Sk,N- Let <fi and tp be Jacobi forms in Sk,N- Then one has 

Proof. Indeed, if x and y are the coordinate vectors of </> respectively tp with respect 
to the basis Ph(i) we have ((f>\ip) = x t Thy- On the other hand, we deduce from 
{(p\Ph(i)) = C^(h(i)) the identity T h x = ^ h {4>) and similarly T h y = j h {ip)- Using 
T/j = the lemma is now immediate. □ 

3.2. Fourier coefficients of Jacobi Poincare series. According to Theorem 1 
the computation of the iV-th coefficient in the series Dp,a is reduced to the numeri- 
cal computation of the Fourier coefficients of the Jacobi Poincare series Pa,t in £fc,iv 
where A < and r 2 = A mod AN. Explicit formulas for these coefficients were 
originally worked out in [8, Prop, in Sect. 2]. However, we shall use the formulas 
from [3, Thm. 1.4], which are stated in terms of vector valued modular forms instead 
of Jacobi forms. The formulas given in [3, Thm. 1.4] are more useful for us mainly 
because they involve in an explicit manner the Weil representation p associated to 
the lattice L := (Z, (x, y) ^ 2Nxy). Explicit formulas for the matrix coefficients of 
this representation have been given in [14] and will be given for arbitrary lattices 
and in a slightly more explicit form in [18]. This makes it possible to evaluate the 
coefficients of the Poincare series without having to deal with too many iterated 
sums 2 (more precisely, Gauss sums inside Kloosterman sums) as provided by the 
formulas in [8]. 

To be more specific, for our computations we shall use the following formulas: 



A generalization and implementation of this method for arbitrary Weil representations will 
be published elsewhere [18]. 
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Theorem 2. We have 



PA,r (z, t) = ajv,fe,A,r Y] PA.r ( A', r') q ' ™ ( r 



r',A'eZ,A<0 
r /2 =A' mod 4AT 



where the Fourier coefficients are given by 

PA,r(A',r') = 6a,A' (o~ r ,r' + S-r,r') 



+2tt 



k 

A' 2 " 



_ ff:(r,A,r',AV,_|(^#Al)- 



i/ere J„ (x) is i/ie standard J-Bessel function and 

mir^.^^t e «(-^)- 

d mod c 
gcd(c,rf) — 1 

where p is the Weil representation o/SL(2, Z) associated to the lattice (Z, (x,y) n- 
2Nxy), and where, for an integer r, we use r for the element ^ +Z m t^-Z/Z (see 
Section 3.3 for an explanation of these terms, and see Theorem 3 for an explicit 
expression for the matrix elements of p ), and where, for each d relatively prime to 
c, we choose an element (" d) * n SL(2,Z). Finally, one has 

(7r|A|/iV) fe - 3 / 2 

CtN,k,A,r 



8VNT(k- 3/2) 

Remark. Note that the Fourier coefficients of the Poincare series are real as follows, 
e.g., from [3, Lemma 1.13 and 1.15]. 

Proof of Theorem 2. According to [4, Theorem 5.1] we have a one to one corre- 
spondence between Jacobi forms and vector-valued modular forms. Namely, if <fi is 
a Jacobi form in J^jv with Fourier coefficients C<p(A,r), then the application 

4>^H{4>):= ^2 C 4>( A >P) e p/2N+z 

p mod N 

defines a bijection H : J k ,N — > M k _i L . Here, for a residue class x in ^Z/Z, we 
use e x for the element in the group ring Cf^Z/Z] which is 1 at x and otherwise. 
Moreover, following the notations in [3, Def. 1.2, p. 18], we use M k _i L for the space 
of vector valued modular forms / : H — > C[^Z/Z] which satisfy 

f(Ar)=w(r) 2k - 1 p*(A,w)f(T) 

for all (A, w) € Mp(2,Z) with Mp(2,Z) denoting the usual metaplectic cover of 
SL(2,Z) (see [3, sec. 1.1, p. 15]). The representation p is the Weil representation 
associated to the lattice L = (Z, (x, y) >->• 2Nxy), i.e. the representation of Mp(2, Z) 
on C[27yZ/Z], which in terms of the standard generators for Mp(2,Z) is given by 

p{ ( h I ) , 1) e r/2N+z = e 2- 2 /4A 

(5) „-2iri/8 

P\\i o ) > V e r /2A+z = 2^ e e rV2A r +z 

r' mod 2JV 

Finally, p* is the representation obtained from p by taking the complex conjugates 
of the p(A, w) (where we view p(A, w) as a matrix). 
By [4, Thm. 5.3] we have 

(6) m) = ^=(H(<f>)\H(1>)h 
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where the scalar product of two vector valued modular forms /i and fa in M k _i L 
is defined as [3, cq. (1.17), p. 22]: 

(/i|/ 2 )= I h,j(r)f 2 , j (T)v k -*dudv {fj= h^ r ). 

(Note that [4, Thm. 5.3] contains an error: the factor l/v / 2m in Thm. 5.3 has to 
be replaced by 1/V4m as the proof in [4] actually shows.) In addition, one has to 
use that a fundamental domain for SL(2, Z) J \H x C is given by {(r, xt + y) : r € 
D, x, y > 0, x + y < 1}. Here D is a fundamental domain for SL(2, Z)\H. Note the 
condition x + y < 1, which is due to the action of the scalar matrix — 1 in SL(2, Z), 
which gives (— l)(r, z) = (r, — z). This gives another factor of \, so that the correct 
factor in Thm. 5.3 becomes in fact 1/4-y/m (instead of \j\J1m). 
From (6) and [3, Prop. 1.5, p. 23] it follows that 

H{P&,r) = a N,k,A,rPr/2N+Z,\A\/4Ni 

where P^/ 2 n+z |A|/4JV ^ s tnc P° mcar e series as introduced in [3, Eq. (1.11), p. 19]. 
The Fourier expansion of P^n+z |A|/4JV nas been computed in [3, Thm. 1.4, p. 19]. 
The theorem is now an immediate consequence of the formula loc. cit. □ 



3.3. The Weil representation. The last requirement to make our formulas for 
the Petersson scalar products of Jacobi forms explicit and ready for implementation 
is the determination of the matrix elements of the Weil representation p which is 
defined by the formulas (5). Explicit formulas are given e.g. in [13, Prop. 1.6]. 
However, a direct implementation of these formulas would involve, for A = ("J) 
with c / 0, a sum of length |c|. A trivial estimate shows that with this formula 
the number of exponential evaluations necessary to compute pA.r (A', r') is of the 
order of magnitude M 3 if we truncate the outer sum (5) at |c| = M. The formula 
of Theorem 3 below which we will use is computationally less expensive since it 
requires only one exponential evaluation instead of a sum of |c| such evaluations. 

For stating our formula we need some notations. To each matrix A = ( ° b d j in 
SL(2,Z) we associate an element A* in Mp(2,Z) by setting A* := (A, \J ct + d). 
(Note that A i— )• A* does not define a homomorphism of groups.) Here and in the 
sequel, for a complex number w, we use y/w for its root which is situated in the right 
half plane but not on the negative imaginary axes. We write e n (w) and e„(u>) for 
e 2TTinw anc j e 2TTiw/n^ respectively. For h, h' in ^Z/Z we define the matrix element 
p{A*)h,h> by the equation 

p{A*)e h ,= P( A *)h,h>e h . 

Finally, for a nonzero rational number r = - (gcd(a, b) = 1, a > 0) we define 
signature(r) as the integer modulo 8 such that 

es(signature(r)) = [ — ] e 8 (l — a\ + even(a) ba\), 
W 

where even(a) equals 1 or accordingly as a is even or odd, where we write a = 2 n a\ 
with odd ai, and where denotes the generalized Legendre symbol (which is 
multiplicative in a and in 6, and which satisfies (|) = 1 if b = ±1 mod 8 and 
(|) = -1 if b = ±3 mod 8. 
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Lemma. Let r be a nonzero rational number, a its (positive) denominator, and 
suppose that a is not exactly divisible by 2. Then 



— j= e(rx 2 ) = e 8 ( signaturc(2r)) . 

V5 



x mod a 

where a = a/ gcd(2,a). 

Proof. This identity can be easily verified using the well-known formulas 

'b 



V" » a (bx 2 ) = I (%)iy/a ifa = 3mod4 

zmoda e 2mb ^V^a if a = 2™, n>2. 

for standard Gauss sums [7], and substituting x = 2 n t + a\u in a Gauss sum for 
a general a = 2 n a\ (with odd a\), where t and u run through a complete set of 
representatives modulo oi and 2", respectively. □ 

Theorem 3. For A = ( a c h d ) in SL(2,Z) and integers x, x' , the matrix ele- 
ment p{ A*) h . h , (h = x/2N + Z, ft' = x'/2N + Z) of the representation p o/Mp(2,Z) 
defined by (5) is given by the following formula: 



(£) Va if a = 1 mod 4 



e 



p(^*)fc,fc' = x(A)y g ° d ^' ^ e 4JV (a^ 2 + 26ca;y + cdy 2 + 2z c (bx + dy) + abz 2 c ) 

if there exists a y such that x' = ax + cy + z c mod 2N ( and then the right hand 
of this formula does not depend on the choice of y), and p{A*)h,h' = otherwise. 
Here z c — N if c/2N is a 2-adic unit, and z c — otherwise. Finally, 



X(A) = e 8 (-sign(c)) 



e 8 ( signature(aiV/2c)) if c/2N is a 2-adic unit 
e 8 ( signature(2a./V/c)) otherwise. 



if c ^ and x(^) = l/\/sign(d) otherwise. 

Proof. The claimed formula is proved in [14, Sec. 1.4] (see formulas (26), (26)' and 
(27) on p. 37 loc.cit.), and in a slightly different form in [17, Prop. 4.1]. Note that 
the formulas in [14] refer to a right action a) i-> $\ i N a of Mp(2, Z) on a certain 
space ThN spanned by theta functions "&n, p (p mod 27V). The precise definitions 
of these objects is not important here. For applying the formulas in [14] one has to 
identify C[^Z/Z] with the dual of the space Tftjv via eh($N.p) = 1 if ft = p/2N+Z 
and eh{&N, P ) = otherwise. It is then easily verified (using [14, Satz 0.1, p. 10]) 
that (p(a) eh) ($) = e/ l (t?|i Jv a) for all ft in -^ZjZ and all i9 in TftM- Inserting 
for ??|i !Ar a: the formulas (26), (26)' and (27) of [14, Sec. 1.4] proves the claimed 
formulas apart from the fact that the quantity x(A) is given in [14] in the form 

'l/yiig^d) c = 

X(A) = Vgcd( 1 2W ,c)c, mod 2c e 4c (aiV^) if c/2AT is a 2-adic unit 
k £, mod c e c («iV, 2 ) otherwise. 

(Note that formula (27) contains a typo stating that x( a ) = 1 for c = 0, which, 
however, can be easily corrected by noticing that $N,p\± t N ( ~o -l) = i~ 1 '&N,- P -) 
By the Lemma we can rewrite x(A) m the given form. This proves the theorem. □ 
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3.4. Implementation. All the computations were done using Sage [12]. The rou- 
tines for interval arithmetic are implemented in Sage via the package mpmath. The 
code for computing the Weil representation, the Poincare series and finally the files 
containing the necessary Gram matrices Th can be found on [11]. The numerical 
computation of the Dirichlet series Dp t c(s), which needs in addition to the Gram 
matrices the Fourier coefficients of the eigenforms in question is done via a short 
Sage script (also to be found loc.cit.), which in turn calls the Siegel modular forms 
package [10] mentioned in Section 1. 

3.5. Discussion of error estimates. For the computation of the Gram matrix 

of Theorem 1 we need to compute the Fourier coefficients pA, r (A', r') of Theorem 2. 
Using the standard asymptotics of the J-Bcsscl function J v (x) ~ r(v+i) (f )" + 
for x <C 1 and the fact that \p (M) h h , | < 1 it follows that in order to obtain an 
error of order e = 10~ T for the computation of a coefficient we need to truncate 
the sum over c at about c — M (T) with 

T 

M (T) > O(l)10^ 

with a constant (which can be made explicit) depending on fc, A, A' and N. From 
this we see that for small values of k, in particular, this evaluation can be very 
time consuming while for higher weight k we can truncate the sum earlier. At 
this point one should also mention that the dimension of the space Sk,N, i.e. the 
number of necessary Poincare series to consider, increases linearly in k and in N 
(see [16, Thm. 1] for the dimension formula), so that the number of coefficients 
PA,r (A',r') to compute grows quadratically in k and in N. In Table 2 we list the 
dimensions of the spaces of Jacobi cusp forms S k ,N which fall within the range of 
our computations. 

We need to make sure that the set of Poincare series {■?*/,(») for a given 

sequence of indices h = {(A i ,r i )} i (n — dimSj^jv) form a basis of Sk,N- For this it 
suffices to verify that the associated Gram matrix Th (4) is invertible, i.e. that its 
determinant is non-zero. Let e = 10~ T (T 1) and suppose that we use a suffi- 
ciently high working precision. Let N (Th) — ( a ij) be a numerical approximation 
to Th with an error less than e. By this we mean that = N(Cp h[i) (h (j))) and 
maxo<i,j<„ \ctij - C Ph(i) (h(j))\ < e. It follows that -e < - C Ph(i) (h(j)) < e, 
or equivalently Cp h(i) (h(j)) G (a^ — e, a^- + e) =: hj for each Let I (Th) be 
the interval matrix with intervals Iij as entries. The matrix / (Th) is then also an 
approximation to Th in the sense that I\ e I (Th)- By elementary properties of 
interval arithmetic (see e.g. [1]) one can see that det Th G det I (Th) and it is there- 
fore enough to verify that ^ det I (Th)- In practice, we start with one function 
Ph(i) an d the corresponding one-by-one gram matrix, T h X \ We then choose a sec- 
ond function, Ph(2), an d calculate the determinant (as an interval) of the interval 
approximation to the partial Gram matrix of these two functions. If this interval 
does not contain zero we add the function Ph(2) to the set and have a non-singular 
two- by- two matrix T y h . This procedure then continues until we have a set of n 
linearly independent functions and a non-singular matrix IV 

It is not difficult to see that if the absolute error in (the infinity norm of) of 
our approximation N(Th) of Th is e then the absolute error of our numerical ap- 
proximation N({Fn\Gn)) of (Fn\Gn) accordingly is also bounded by a constant 
(depending on Fn and Gjv) times e. To see this, first observe that in all our ap- 
plications the functions Fn and Gm are known exactly (if defined over a number 
field other than Q we can approximate them to arbitrarily small precision). For 
any tp,ip e S k ,N let (<p\ip)N(r h ) : = lh (</>)* N (T h ) -fh (V>) denote the inner-product 
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with respect to the matrix N(Th). Then 

(7) <\\Th-th\u\ih{<p)hh h mw 

Thus the absolute error depends on the Z 2 -norms of the vectors jh (ip) and 7^ (tp). 
Since the numbers N{(Fn\Gn}) tend to grow rapidly with TV it makes more sense 
to consider, as a measure of accuracy, the number of correct digits, i.e. the relative 
error, instead of the absolute error. Unfortunately we are not able to estimate 
(</?!?/>) from below and thus can not give a corresponding "predictive" bound for the 
relative error in N((Fn\Gn)). However, it is clear that if 

\N((F N \G N )) - (F N \G N )\ < e ■ \\ lh (F N ) \\ 2 \\ lh (G N ) || 2 < e ■ K 

and \e ■ K\ < \N({Fn\Gn})\ then the relative error can be estimated by 

\N((F N \G N ))-{F N \G N )\ e_K_ 

\(F N \G N )\ \\N((F N \G N ))\-e-K\ 

and we are therefore able to give a precise a posteriori bound for the relative error. 
Provided that N {(Fn\Gn)) is nonzero (which turns out to be the case in all cases 
we considered so far) we are always able to go back and decrease the initial e until 
we have achieved any desired bound. Since rigorous, or even validated, numerics 
never give insurance against e.g. errors in the actual implementation we also felt 
compelled to apply our algorithms to a case where we know the results explicitly. 

Recall that if G is in the Maass spezialschaar and F is an cigenform, then 
Df,g(s) — (F\\Gi)L(F, s)L(F, s). To verify that the Gram data is correct as well 
as that the relative error behaves well, we compare the computed values for the 
coefficients of Df.g( s ) with those of L(F, s) (see Tables 3 and 5). The coefficients 
of L(F,s) can be computed exactly using [10], It seems that in most cases the 
relative error can be estimated from above by a constant (not growing very fast 
with respect to N) times — f^Hcx,. That is, the quotient of N{(Fjv\Gn)) divided 
by \\lh {Fn) II2II7/1 {Gn) H2 does not seem to tend to zero too fast with N. 

Another way we verify that the Gram data is correct is we check, as in [6] , that 
the series Dfg for F and G not in spezialschaar, satisfies the functional equation 
shown in [9]. 

3.6. Runtimes. To give a flavor of the amount of time needed to run an actual 
computation of the Gram matrices we can consider weight k = 30 where we spent 
slightly less than 4h (14254s) on a 2666MHz processor to compute the Gram ma- 
trices for 1 < N < 20 with an error bound for the Poincare series of 10 -50 . This 
involved computing a total of 16340 Fourier coefficients. We note that the time 
for a single Poincare series computation to a given precision is essentially bounded 
by a constant independent of k and N. Thus the runtime only increases with the 
number of coefficients required. The dependence of the runtime on the precision 
requested is slightly more complicated. For the proof of the Main Theorem in the 
case k — 30 we needed a precision 10 -90 and the corresponding Poincare series 
computations (with a working precision of 631 digits) took more than two weeks 
on the same system as mentioned above. 

4. Proof of the Main Theorem and its Supplement 

Fix a weight k (between 20 and 30), and let F\,. . . , F4 be a basis of common 
Hecke eigenforms of the full space of Siegel modular forms of weight k. We 
assume that F\,. . . , Fh are the interesting ones, i.e. that they span the space S\. 
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We want to show that the series Dp i ,p j (s) and L(Fi,s) (1 < i < j < h, 1 < I < d) 
are linearly independent. For this we set 



M := 



( (-Fi,i| F i,i>/7i,i 

<-Pl,l |-P 2 ,1>/T1,2 



<F e ,i|F e ,i>/ 7cje 
A fi,i 



(Fi, 
(Fi, 



\Fi,, 
\F 2 ,, 



>/7i,i\ 

>/7l,2 



<F e ,„|F e ,„}/ 7e 



V 



Here n 



e(e+l) 



d, we use -F^jv for the A^-th Fourier- Jacobi coefBcient of Fi, 
and we use jij — (Fi,i\Fj t i). The Xp. , are the eigenvalues of the Fj with respect 
to T{1) (see (1)). We want to show detM ^ 0. The eigenvalues Xp jl are totally 
real numbers and can be computed exactly [15]. Assume we approximate them 
numerically by rational numbers N(Xp j ; ) such that 



IA, 



N(X F . t )\<e. 



Furthermore, using the method described in Section 3 we compute rational approx- 
imations N(((Fi,i\Fj,i) /-fij) such that 



(8) 



\(F i , l \F Jtl )-N{(F itl \F j , l ))\/\>y i j\<e. 



(The numerical calculations show that indeed jij ^ for all If we denote 

by N(M) the matrix obtained from M by replacing the coefficients m p ^ q by their 
respective rational approximations N(m Pyq ) we conclude that \\M — N(M)\\ OQ < e. 
As in section 3.5 we construct an interval matrix approximation to M, I(M), with 
intervals 1 VA = (N (m p ^ q ) — e, N {m p ^ q ) + e) as entries. Since M e I(M) we are 
once again able to show that the determinant of M is non-zero by verifying that 
the determinant D (which is an interval) of I(M) does not contain zero. If we 
compute the gram matrix T/j to the precision e\ we obtain an upper bound e 2 for 
the relative error e in (8) by 



( 2 = £i 



max; 



hh(Fi, 



mm; j 1 7i j I 

It turns out that by choosing ei as in Table 6 we obtained an e 2 such that the 
determinant of the corresponding interval matrix I(M) did not contain zero. This 
proves the Theorem. 

For the actual computed values of upper and lower bounds of ||7h(^i,j)ll2 and 
|7i,j|, respectively, see Table 5. The values of e 2 obtained from these bounds, as 
well as the determinant D = (D m — Dg, D m + Ds) are given in Table 6. That zero 
is not contained in either of these intervals is easily checked. It should also be 
remarked that for all our computations for a given fc, we used a working precision 
set high enough to represent D with an absolute error comparable to e 2 . As an 
example we computed the Poincare series for k = 30 using 631 digits precision. 

The supplement was proved in a similar manner, using interval arithmetic and 
the error bounds (7). 



5. Concluding Remarks 

We computed Df,g{ s ) numerically following the method of Section 3 for those 
cusp forms F and G listed by name in Table 1. The Fourier coefficients of these 
Siegel cusp forms and the coefficients of their spinor L-functions can be computed 
exactly using the package [10] which implements the methods of computation of [15] . 
As additional control for our calculations, we verified that our numerical approxi- 
mations of Dp.a, for F as in Table 1 and for a Maass spezial form G, equals indeed 
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(within our error bounds) the spinor L-function of F, as predicted by the second 
theorem in [9] (see also Tables 3 and 5). 

Andrianov asks [2, Concl. 2, p. 114] if all Dirichlet series satisfying the functional 
equation (2) are linear combinations of the spinor L-scries L(H, s) where H runs 
through the eigenforms in the full space Mk of Siegel modular forms of weight k 
on the full Siegel modular group. We can answer this question in the negative. 
Namely, we computed all the Hecke eigenforms in weights 20, 22, 24, 26, 28 and s30 
and their corresponding L-series; in particular we computed the L-series for the 
Eisenstein series, Klingen-Eisenstein series, Maass lifts, and interesting forms that 
are Hecke eigenforms. For all of the cases listed in Table 1, Df^g(s) was not a linear 
combination of spinor L-series. 

This raises the question of the nature of the series Dfq for interesting forms F 
and G. It is natural to ask first of all whether they are L-functions at all in the 
sense of Selberg, i.e. whether they belong to the Selberg class. In fact, they do not, 
since they do not possess an Euler product expansion (see the supplement to the 
main theorem and Table 4). If we set 

(9) D F , G (s) = (F^G^Dp^s + k - |). 

(the scalar product (F\\Gi) is different from in all cases covered by Table 1) the 
functional equation (2) becomes 

D F , G (s) := (2n)- 2s - 2k + 3 T( S + k ^)T(s + \)D F , G {F, s) = D^ G (1 s). 

Accordingly to [9, Theorem 1] the function Df,g(s) has a pole at s — |. This 
shows again that ^^^(s) is not in the Selberg class. 
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6. Tables and Figures 

Table 1. Dimensions of subspaces of M^. (KE = Klingen- 
Eiscnstein series, MC = Maass spezialschaar cusp forms, Eigen- 
forms = names of Galois orbits of interesting eigenforms following 
the naming conventions in [15]) 



k 


Sk 


KE 


MC 


sl 


Eigenforms 


20 


5 


2 


2 


1 


"f~20 


22 


6 


2 


3 


1 


T 22 


24 


8 


3 


3 


2 


i24a; "^246 


26 


7 


2 


3 


2 


T26a, T 2 6b 


28 


10 


3 


4 


3 


T 28 


30 


11 


3 


4 


4 





Table 2. Dimensions of the spaces Sk,N for even weights 20 < 
k < 30 and indices 1 < N < 20. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


20 


2 


4 


5 


6 


8 


9 


11 


12 


13 


15 


17 


17 


20 


21 


22 


23 


26 


26 


29 


29 


22 


3 


4 


6 


7 


9 


10 


13 


13 


15 


17 


19 


19 


23 


23 


25 


26 


29 


29 


33 


32 


24 


3 


5 


7 


8 


10 


12 


14 


15 


17 


19 


21 


22 


25 


26 


28 


29 


32 


33 


36 


36 


26 


3 


5 


7 


8 


11 


12 


15 


16 


18 


20 


23 


23 


27 


28 


30 


31 


35 


35 


39 


39 


28 


4 


6 


8 


10 


12 


14 


17 


18 


20 


23 


25 


26 


30 


31 


33 


35 


38 


39 


43 


43 


30 


4 


6 


9 


10 


13 


15 


18 


19 


22 


24 


27 


28 


32 


33 


36 


37 


41 


42 


46 


46 



Table 3. For the Maass lift G of the first Fourier- Jacobi coeffi- 
cient <j) °f Y24a, we compute numerically the coefficients a(N) of 
-Dr24a,G( s )/(0|0) by the method in Section 3 and we compute ex- 
actly the coefficients b(N) of L(T24a, s) (which are integers). This 
is an additional check that our method of computation gives cor- 
rect results since a theorem in [9] tells us that, in this case, these 
two Dirichlet series coincide. 



N 


a(N) 


b(N) 


1 


1.0000000000000000000000000000 




1 


2 


-5.5603200000000000000000000000 


10 6 


-5560320 


3 


-5.3017924680000000000000000000 


10 10 


-53017924680 


4 


4.3592282275840000000000000000 


10 12 


4359228227584 


5 


-3.3324163624500000000000000000 


10 13 


-33324163624500 


6 


2.9479662695669760000000000000 


10 17 


294796626956697600 


7 


8.9548405531223548000000000000 


10 18 


8954840553122354800 


8 


-1.7002266475219648512000000000 


10 20 


-170022664752196485120 


9 


3.3581886607436193369000000000 


10 19 


33581886607436193369 


10 


1.8529301348457984000000000000 


10 20 


185293013484579840000 
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Table 4. The first 18 coefficients of Dx24x,T24y(s) (see(9)), where 
x,y e {a,b}, where T24a, T246 denote the two interesting eigen- 
forms in weight 24. Note that Dr24x.T24y(s) satisfies a functional 
equation with respect to s H 1 - s. More data and this data with 
more accuracy are available at [11]. 





-DT24a,T24a(s) 


-DT24a, T24b{s) 


-C>T24b,T24&(s) 


1 


1.0000000000... 


1.0000000000... 


1 


0000000000 . . . 


2 


1.2562996887... 


-1.6301379426... 


1 


1807313893... 


3 


1.7810603106... 


2.0303021423... 





8677710052 . . . 


4 


2.0741433142 . . . 


-3.4104218951... 


2 


8372669032 . . . 


5 


2.8899783797... 


0.0053403337 . . . 


2 


3422724376 . . . 


6 


2.6795441187... 


3.4130123501 . . . 


2 


5830326424 . . . 


7 


3.6002676445 . . . 


1.2502722224... 


2 


2422068960 . . . 


8 


3.0295878336 . . . 


-0.2164387720... 


3 


5760646077 . . . 


9 


3.9970248055 . . . 


-1.9659335003... 


3 


2031268531 . . . 


10 


2.9874427387 . . . 


-2.2751217843... 


3 


5506087520 . . . 


11 


3.5420885329 . . . 


-0.5800680650... 


2 


6657794870 . . . 


12 


2.9887639258 . . . 


2.3742550783 . . . 


3 


3541225360 . . . 


13 


4.7252631756... 


-1.4101441627... 


3 


0042305971 . . . 


14 


3.6881022526... 


2.8085205331 . . . 


2 


5295615640 . . . 


15 


3.2033969064... 


-2.1277658831... 


2 


9089500309 . . . 


16 


4.8466409087 . . . 


1.6172735165... 


5 


3475003622 . . . 


17 


5.5339720201 . . . 


-2.0014367544... 


4 


1580587188... 


18 


4.5912959032... 


-2.3762191782... 


4 


6382023273 . . . 



Table 5. Computed upper bounds for the I 2 norms of -fh (Fij) 
and lower bounds for the inner products 7^ = (Fi t i\Fj t i). The last 
column contains the maximum error in the comparison between the 
Hecke eigenvalues Xp^j and the numerical approximations, Af^ , 
computed using the inner product of Fi with a form from the Maass 
spezialschaar. In all cases 1 < i, j < n where n, given in the second 
column, is the dimension of the space of L-series which we want to 
show is linearly independent. 



k 


n 


max(|| 7ft (F u )\\ 7 ) 


min \-fij\ 


max \X FiJ - X FiJ | 


20 


5 


1.2 • 10 U9 


1 • 10 1 


3.2 • 10~ 45 


22 


6 


8.2 • 10 U4 


5 • lO" 11 


9.6 • 10~ 5ti 


24 


11 


7.8 • 10 U9 


1 • io- iu 


1.7 • 10" bb 


26 


10 


5.2 • 10 1U 


1 • 10" 12 


7.5 • 10-' 2 


28 


16 


8.8 • 10 ib 


1 • 10° 


2.2 • 10~ Yi 


30 


21 


2.9 • 10 12 


3 • 10~ 25 


5.6 • lO" 96 
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Table 6. The values used for e\ and the computed values of 62 as 
well as the computed values of the interval determinant of I(M), 
D=(D m -D 5 ,D m + Ds). 



k 


ei 


£2 


D m 


D S 


20 


1 • 10- 4,) 


7.9 • 10~ 26 


-7.05 • 10 054 


1.26 • 10 033 


22 


3 • 10- b0 


1.2 • 10" 23 


-6.63 • 10 usu 


8.87 • 10 ueu 


24 


2 • 10- 5() 


2.8 • 10- 13 


7.38 • IO 209 


1.63 • 10 2Ob 


26 


7 • 10~ bb 


4.0 • 10-' A ' 


-6.80 • 10 iyu 


1.27- 10 iY3 


28 


1 • 10" e5 


4.8 • IO- 33 


-3.21 • 10 399 


6.45 • 10 384 


30 


1 • io- yu 


2.9 • 1Q- 2 ' 


1.53 • 10 bUY 


1.13 • 10 bUY 
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